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ABSTRACT 

Sub-mm detected galaxies located at redshifts z > 1 host a major fraction of the bolometric luminosity 
at high redshifts due to thermal emission from heated dust grains, yet the nature of these objects remains 
a mystery. The major problem in understanding their origin is whether the dust heating mechanism is 
predominantly caused by star formation or active galactic nuclei, and what triggered this activity. We 
address this issue by examining the structures of a sample of 11 sub-mm galaxies imaged with STIS on 
the Hubble Space Telescope. We argue that ~ 61 ± 21% of these sub-mm sources are undergoing an 
active major merger using the CAS quantitative morphological system. We rule out at ~ 5 a confidence 
that these sub-mm galaxies are normal Hubble types at high redshift. This merger fraction appears 
to be higher than for Lyman-break galaxies undergoing mergers at similar redshifts. Using reasonable 
constraints on the stellar masses of Lyman-Break galaxies and these sub-mm sources, we further argue 
that at redshifts z ~ 2 — 3 systems with high stellar masses are more likely than lower mass galaxies to 
be involved in major mergers. 
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1. INTRODUCTION 

Faint sub-mm galaxies were discovered with the first 
generation of multi-element detecting devices working at 
sub-mm wavelengths, most notably the SCUBA array 
(Smail, Ivison & Blain 1997). The nature of these galaxies 
has however remained a mystery, despite a considerable 
amount of observing time spent obtaining ~ 300 detec- 
tions in the field and behind lensing clusters (Blain et al. 
2002). 

The generally accepted working idea is that these sub- 
mm galaxies are distant z > 1 systems which emit in the 
rest-frame far infrared due to thermal emission from dust 
grains heated by star formation and/or active galactic nu- 
clei (AGN). These objects are also thought to be analogs 
of nearby ultra luminous infrared galaxies (ULIRGs), al- 
though even in nearby ULIRGs it is debated whether or 
not the dust is heated by the UV continuum of massive 
stars, or an active nucleus. As sub-mm galaxies are 400 
times more common at z ~ 2 than at z ~ 0, they con- 
stitute a significant fraction of the bolometric luminosity 
density at high redshift (Chapman et al. 2003a), thus un- 
derstanding their origin and relationship to nearby galax- 
ies is of central importance. 

Answering fundamental questions concerning sub-mm 
galaxies by studying them at multiple wavelengths has re- 
mained very difficult due to the low-resolution of SCUBA 
and other sub-mm instruments, making identifications at 
other wavelengths difficult. One approach to this prob- 
lem has been to identify sub-mm galaxies through their 
emission in the radio (Ivison et al. 1998; Chapman et al. 
2002a), utilizing the fact that sub-mm/far-infrared lumi- 
nosities correlate in nearby galaxies (Helou et al. 1985). 
Through these identifications, we are able to determine 
sub-mm source positions to within a sub-arcsec, from 



which optical follow-up can be done readily. A subset of 
these sub-mm galaxies whose positions were located using 
the radio have been imaged with the Hubble Space Tele- 
scope using STIS (see also a companion paper, Chapman 
et al. 2003b). Previous studies have been largely quali- 
tative, and often with misidentifications (e.g., Smail et al. 
1998; Ivison et al. 2001; Chapman et al. 2002b). 

Based on our analysis we conclude that 40 - 80% of 
these sub-mm sources are consistent with undergoing a 
major merger, and statistically rule out at 5 a confidence 
that these systems are normal galaxies at high redshift. 
We further demonstrate, through comparisons with artifi- 
cially redshifted nearby ULIRGs and normal (non-ULIRG) 
galaxies, that high redshift sub-mm galaxies appear qual- 
itatively similar to ULIRGs, and are potentially forming 
into massive spheroids. 

2. IMAGING AND ANALYSIS METHOD 

Our sample consists of 11 sub-mm sources selected in 
the radio. The full sample selection for these galaxies is 
described in Chapman et al. (2003b). Although there are 
some selection issues based on using radio identifications, 
these objects were not selected based on optical properties, 
and they span a wide range of optical magnitudes (Chap- 
man et al. 2003b). Each of the sub-mm sources (Figure 1) 
we study was imaging in the 50CD clear filter in one to 
three orbits, in two separate exposures per orbit, result- 
ing in a total exposure time of 1-3 ks (see Chapman et al. 
2003b for further details). 

We analyzed these images using the CAS (concentra- 
tion, asymmetry and dumpiness) morphological system 
(Conselice 2003). In the CAS system, the asymmetry in- 
dex (A) is used to determine whether or not a galaxy is 
involved in a major merger (Conselice et al. 2000a, b; Con- 
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selice 2003; Conselice et al. 2003). The value of A is cal- 
culated by rotating a galaxy through 180° and subtracting 
this rotated galaxy from the original and comparing the 
absolute value of the residuals of this subtraction to the 
original galaxy flux (Conselice et al. 2000a). Before car- 
rying out the CAS analysis on the sub-mm galaxy images 
we removed any galaxies and stars close to the sub-mm 
sources, so as to not contaminate the measurements. We 
placed our initial guess for the center on the brightest por- 
tion of the sub-mm galaxy, and then ran the CAS program 
to determine asymmetries and light concentrations. 

To understand the systematics due to a lowered resolu- 
tion and signal to noise we artificially redshifted 50 nearby 
ULIRGs (from Farrah et al. 2001 and Conselice 2003) 
and 82 normal Hubble types (Frei et al. 1996, Conselice 
2003) to how they would appear in our STIS images at 
z ~ 2 — 3, the spectroscopic redshift range of these sources 
(Chapman et al. in prep). These simulations are done by 
reducing the resolution and surface brightness of a galaxy 
as to how it would appear at high redshifts. The amount 
of noise matching that expected in the STIS observations 
are then added in. 

Once we create these new images, based on these simula- 
tions, we remeasure the CAS parameters of the simulated 
galaxies, in the same manner as on the original sub-mm 
sources. Note that we are taking a very conservative ap- 
proach by including these two types of nearby galaxies to 
account for redshift effects. The apparent morphology of 
the sub-mm sources is very peculiar and therefore the like- 
wise peculiar nearby ULIRGs are a better population for 
making this correction. Normal galaxies also do not be- 
come significantly more irregular in the rest-frame ultra- 
violet (Windhorst et al. 2002). We use the nearby normal 
sample as an extreme lower limit to what these corrections 
could be. 

3. RESULTS 
3.1. Evidence For a Merger Origin 

Concentration-asymmetry diagrams for our sample of 
sub-mm galaxies are shown in Figure 2. The asymmetry 
and concentration values for these galaxies have been cor- 
rected for redshift effects assuming they have morphologies 
intrinsically similar to ULIRGs (Figure 2a, c), and nearby 
normal galaxies (Figure 2b, d) (Conselice et al. 2000a). 

We also plot on Figure 2a,b the average values and lcr 
variations of measured concentrations and asymmetries for 
nearby galaxy populations, including ULIRGs, as observed 
in the rest-frame optical (Conselice 2003). Although we 
view these sub-mm galaxies in the near to mid-UV, the 
morphological appearance of ULIRGs and other star form- 
ing galaxies does not significantly differ between optical 
and near-UV wavelengths (e.g., Conselice et al. 2000c; 
Surace & Sanders 2000; Windhorst et al. 2002). Fig- 
ure 2c, d shows the location of the sub-mm galaxies in ref- 
erence to the actual C — A values found for nearby nor- 
mal galaxies, including separately labeled ellipticals and 
ULIRGs. 

Figure 2 can be used to determine which population the 
sub-mm galaxies are most similar to morphologically. As 
can be seen in Figure 2, the dominate morphological fea- 
ture of the sub-mm galaxies is their high asymmetries. The 
distribution of the sub-mm sources in C — A space is also 



most similar to the ULIRGs, and does not overlap much 
with any normal galaxy type. Using the major merger 
criteria calibrated in Conselice (2003), based on nearby 
ULIRGs, we find that a galaxy is likely a major merger if 
A > ylmorgcr = 0.35. Using the ULIRG z ~ 2 correction, 
we find that the merger fraction for the sub-mm sources is 
0.82, which is the same fraction found when correcting by 
the z ~ 3 simulations. The normal galaxy correction still 
reveals a large merger fraction of ~ 0.4, which is certainly 
a lower limit to the actual fraction of galaxies involved in 
mergers. This lower limit is as high as the largest merger 
fractions found for any Lyman-Break galaxy population 
(Conselice et al. 2003)' (§3.2). 

There are naturally uncertainties in this merger frac- 
tion calculation. First, the sample size is small, only 11 
galaxies. Second, it is based on a correction which has nec- 
essarily some uncertainty in it, since the sub-mm galaxies 
may be a mix of normal and mcrgcr/ULIRG systems. To 
circumvent this problem we performed a series of Monte 
Carlo simulations using our normal and ULIRG galaxy 
samples after they were simulated to z ~ 2. We take at 
random 11 of these objects and compute their asymmetry 
distributions, and compare this to the asymmetry distri- 
butions of the sub-mm sources. Doing this, we find at a 
5 a confidence that normal galaxies at z ~ 2 cannot repro- 
duce the asymmetry distribution of the sub-mm sources. 
The simulated ULIRGs are also slightly less asymmetric 
than the sub-mm sources, but have asymmetry values dis- 
tributions consistent to within 1.9 a. 

3.2. Sub-mm and Lyman-break Galaxies 

The implied merger fraction of these sub-mm galaxies 
appears to be slightly larger than what is found for UV 
bright galaxies in the Hubble Deep Field (HDF) (Conselice 
et al. 2003) at similar redshifts. Lyman-break galaxies at 
z ~ 2.5 found in the HDF all have implied merger fractions 
lower than our derived sub-mm values (Conselice et al. 
2003). For example UV bright galaxies at 2 < z < 3 with 
magnitudes Mb < —20 and stellar masses > 10 9 5 M have 
a merger fraction ~ 0.18 (Figure 3). The brightest and 
most massive galaxies seen in the HDF, with Mb < —21 
and M* > 10 10 M Q , have merger fractions ~ 0.4 — 0.5, 
which is lower than the implied merger fractions of our 
sub-mm sample. From this it appears that the sub-mm 
galaxies are actively undergoing mergers in a greater abun- 
dance than Lyman-break galaxies. This may be the result 
of higher mass systems undergoing more mergers at higher 
redshifts (Figure 3). Sub-mm galaxies as a population ap- 
pear to be dominated by systems involved in major merg- 
ers, while the Lyman-break galaxies are more likely to be 
in various phases of evolution, and perhaps includes star- 
bursts not induced by major mergers. On Figure 3 we 
plot the merger fraction of the sub-mm sources at a stellar 
mass of 10 11 M Q . This estimate is likely roughly correct 
as sub-mm galaxies have dynamical and gas masses > 10 11 
Mq (e.g., Frayer et al. 1998). 

4. COMPARISON TO Z~ GALAXIES: ULIRGS AND 
SPHEROIDS 

Although we are in the regime of small number statis- 
tics, we can argue that the sub-mm galaxies are similar 
to nearby ULIRGs in terms of their structural properties, 
in addition to their already well established similarities 
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in producing rest-frame far infrared light (e.g., Dey et al. 
1999). 

We have already argued this through the high asym- 
metries and implied high merger fraction for the sub- 
mm galaxies, and the low probability that these sub-mm 
sources are morphological similar to normal galaxies. This 
can be shown qualitatively as well. Figure 4 shows nearby 
ULIRGs at z <~ 0.1, and the same galaxy after it has 
been simulated to how it would appear at z <~ 3. These 
are the same simulations discussed in §2 to determine the 
correction to the concentration and asymmetry values for 
the sub-mm galaxies. These simulated ULIRGs appear 
very similar to the sub-mm galaxies themselves (Figure 1), 
although the separation between components in nearby 
ULIRGs and the sub-mm sources differ (Chapman et al. 
2003b). Nearby normal galaxies, such as disks and ellip- 
ticals do not appear in a similar manner when rcdshiftcd 
out to these same redshifts. 

We can further test how similar the ULIRGs and sub- 
mm galaxies are in terms of their morphologies by per- 
forming a K-S test on their asymmetries and concentra- 
tions. The probability that the sub-mm source asymme- 
tries are taken from the ULIRG population asymmetries 
are 21% after correcting by the ULIRG simulation results, 
and 2% when correcting by the normal galaxy simulations. 
The probability of association with the normal galaxies is 
< 0.001% when using the ULIRG correction and 1% when 
using the normal galaxy simulation correction. The con- 
centration index probabilities of association with ULIRGs 
are 96% for the ULIRG correction and 78% for the nor- 
mal galaxy correction. The corresponding probabilities 
for a normal galaxy association are 12% and 32%, respec- 
tively. Although these K-S tests are not conclusive they 
again rule out that the sub-mm asymmetry distribution is 
similar to the asymmetry distribution of normal galaxies. 

If sub-mm sources are analogs of nearby ULIRGs then 



they are perhaps forming into modern ellipticals. Ellipti- 
cals can be uniquely identified in the CAS system by their 
high light concentrations and low asymmetries. The only 
galaxies that have light concentrations as high as ellipticals 
are ULIRGs (Bershady et al. 2000; Conselice 2003), sug- 
gesting a causal connection. Based on Figure 2 it can also 
be seen that some of the sub-mm sources have light con- 
centrations as high as ellipticals, suggesting that at least 
some of these systems are forming into spheroids. We can 
conclusively rule out that all of these systems are spirals, 
as 4 (36%) have light concentrations too high to be forming 
into disks. 

5. DISCUSSION AND CONCLUSIONS 

Sub-mm galaxies and ULIRGs have comparable 
amounts of rest-frame far infrared luminosities, and large 
masses, suggesting that sub-mm sources result from a sim- 
ilar merger origin. However, direct quantitative evidence 
has been lacking. We argue in this paper that a signifi- 
cant fraction of sub-mm sources are galaxies actively en- 
gaged in major mergers by using the CAS morphological 
system (Conselice 2003). On-going major mergers arc de- 
fined as systems with an asymmetry higher than the lo- 
cally calibrated limit of ^merger = 0.35. Using this limit, 
and correcting for redshift effects, we find that 40-80% 
of sub-mm galaxies are consistent with undergoing major 
mergers. These results have important implications for the 
nature of massive galaxy formation. Sub-mm sources con- 
stitute a significant galaxy population at high redshift, un- 
dergoing mergers, that likely later evolve into the massive 
galaxies seen in the nearby universe. These systems are 
also involved in large amounts of star formation that may 
be induced by the major merger, demonstrating that at 
least for the most massive galaxies, their stars are formed 
through merger induced starbursts rather than collapses. 
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Fig. 1. — Montage images of the 12 sub-mm galaxies imaged with STIS on the Hubble Space Telescope. In this paper, we examine the 
quantitative morphological properties of all but one of these systems whose structure is too faint for an analysis. 
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Fig. 2. — Concentration (C)-asymmctry (A) diagram for the sub-mm galaxies studied in this paper. The triangles show the locations of the 
sub-mm sources on these diagrams after (a,c) correcting their measured values using z ~ ULIRGs simulated to z = 2 and (b,d) correcting 
values using z ~ normal galaxies simulated to z ~ 2. The various crosses in (a,b) show where average nearby galaxy populations fall in this 
space and their 1 a variations. Plotted in (c,d) are the individual C — A values for nearby ULIRGs, and normal galaxies, with the ellipticals 
circles. The horizontal dashed line is the asymmetry limit for major mergers, such that a galaxy is likely undergoing a merger if A > j4 mer ger 
= 0.35. 
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Fig. 3. — Fraction of galaxies between 2 < z < 4 consistent with undergoing a major merger, (f), as a function of stellar mass lower limit. 
The stellar masses of the Lyman-break galaxies (LBGs) were obtained through fits to spectral energy distributions acquired by Papovich ct 
al. (2001) and whose merger properties are discussed in detail in Conselice et al. (2003). The stellar mass lower limit of 10 11 Mq for the 
sub-mm galaxies is estimated based on the large gas and dynamical mass measurements of these galaxies. 



7 




Fig. 4. — Nearby ULIRGs at z ~ 0.1, and images of the same galaxies after they have been simulated to z ~ 3, and viewed under the same 
conditions as the sub-mm galaxies. 



